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We propose a new method to determine the coupling of the Higgs boson to charm quarks, via 
Higgs production in association with a charm-tagged jet: pp —>■ he. As a first estimate, we find that 
at the LHC with 3000 fb“^ it should be possible to derive a constraint of order one, relative to the 
SM value of the charm Yukawa coupling. As a byproduct of this analysis, we present an estimate of 
the exclusive pp —>■ hD^*^ electroweak cross section. Within the SM, the latter turns out to be not 
accessible at the LHC even in the high-luminosity phase. 
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I. INTRODUCTION 

While the Yukawa couplings of the heavy third- 
generation fermions to the Higgs boson can be measured 
at the LHG with a 0(10%) accuracy, see e.g. Ref. [T], 
constraining the diagonal Yukawa couplings of the second 
(first) generation quarks at a level close to the Standard 
Model (SM) expectation is very challenging. An interest¬ 
ing possibility, especially for the second generation, is try¬ 
ing to indirectly access these couplings via the radiative 
decays h — >• M j{Z) [2H5], where is a quarkonium 
stateQ As pointed out in Ref. [8], the exclusive h ^ MV 
decays (V = 7, Z, W) may indeed be accessible at the SM 
level at the LHC and represent a precious source of infor¬ 
mation on physics beyond the SM. In the specific case of 
the charm Yukawa coupling (1/), h should be possible to 
obtain bounds 2-3 times larger than the SM value in the 
high-luminosity (HL) phase of the LHC [5]. These con¬ 
straints are driven mainly by the direct search for h ^ cc 
and, to a smaller extent, also by the indirect sensitivity 
via h —>■ J/dtq. 

In this paper, we propose a new method to measure 
Yc by means of Higgs production in association with a 
charm-tagged jet. A particular advantage of this method, 
compared to the search for h —>■ cc, lies in the fact that we 
probe Yc in production -via the interaction with a charm 
quark from the abundant gc initial state- allowing to 
reconstruct the Higgs from its clean decay modes (/i-J- 
77 or /i —)• WW). This procedure strongly reduces the 
problem of the non-Higgs background, compared to h ^ 
cc. Moreover, requiring a single c-tagged jet in the final 
state allows to adopt high-purity (and low-efficiency) c- 
tag algorithms in order to reduce background (mainly 
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from 6-quark jets), compared to the case of two c-tagged 
jets (as in 6, —)■ cc). 

Compared to the indirect sensitivity to Yc in h ^ 
J/fky, our new method has the advantage of being sensi¬ 
tive to Yc at the tree-level and being based on a process 
that, after charm- and Higgs-tagging efficiencies, yields 
0(1000) signal events at the HL-LHC. For comparison, 
we recall that B{h —)• J/A’-f —)• ^ 10“^, corre¬ 

sponding to 0(10) signal events in pp collisions at 14 TeV 
with 3000 fb~^. The main limiting factor of our approach 
is the theoretical uncertainty on a{pp —?> he), as a func¬ 
tion of Yc. This error could be reduced in the future 
by means of higher-order QCD calculations of the ratio 
o{pp —)■ hc)/a{pp —>■ hb) as a function of 1/ and Y),. 


In principle, the production of the Higgs boson in asso¬ 
ciation with a charm jet (or a charm hadron) can also pro¬ 
ceed via electroweak interactions, with the charm being 
produced by a real or virtual W boson. To complement 
this analysis, and previous studies of exclusive hadronic 
Higgs decays [2[31 ElE], we present here the first estimate 
of the electroweak production of the Higgs boson in asso¬ 
ciation with a single D or D* meson {qq hD^*'>). These 
processes are insensitive to the charm Yukawa coupling 
and could have represented a potential background for 
the extraction of Y).. We have analyzed them in generic 
extensions of the SM, along the lines of Ref. [S]. We 
find that, within the SM, the exclusive electroweak pro¬ 
duction should not be visible at the LHC, even in the 
high-luminosity phase. Moreover, we find that these pro¬ 
cess are not competitive with the corresponding exclusive 
Higgs decays {h —>• MV) as far as generic new physics 
(NP) searches are concerned. 


The paper is organized as follows. In Sect.|TT]we intro¬ 
duce the setup to describe Higgs physics with modified 
Yukawa couplings. The QCD-Yukawa pp —>■ he process 
and the corresponding extraction of Y/ is discussed in 


Sect. HI The exclusive electroweak pp hD^*^ process 


is analyzed in Sect. jlVj The results are summarized in 
the Conclusions. 
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II. SETUP 

Within the SM the couplings of the physical Higgs bo¬ 
son to the fermions are completely determined in terms 
of fermion masses. However, in the presence of NP, a 
misalignment between quark-mass and Yukawa matri¬ 
ces is possible. This can be parametrized in a model- 
independent way by adding the D = 6 operators 

Cj = -^ + ( 1 ) 

to the SM Lagrangian. Here, $ denotes the Higgs 
doublet, parametrized in unitary gauge as $ = 

1 /\/2 ( 0 , h + v)’^, where v corresponds to the vacuum ex¬ 
pectation value ($) = l /-\/2 ( 0 ,u)^, h is the physical 
Higgs field, and qL,UR,dii are the chiral SM-quark dou¬ 
blet and singlets (all quark fields being 3-vectors in flavor 
space). Inserting this decomposition of the Higgs doublet 
into as well as into the SM-like {D = 4) Yukawa terms 

with couplings Y gj^, we obtain the fermion masses and 
Higgs couplings in the flavor basis 

CD-ul\M d-jdR, 

^ ( 2 ) 

u,d '' u,d o 

where the Yukawa matrix Y = ^sm + i Cu,d and the 
mass matrix M = ^(Yg,^ \Cu,d) = ^0^ “ 

Cu,d) are independent parameters. After performing a 
rotation to the mass basis 

= M“,g=diag(m„,m„mO, 

m'^ = Vl = diag{md,ms,mb), 

with Uf = U'l V CKM, we finally arrive at the cou¬ 
plings of the physical quarks to the Higgs boson Y“ = 

= U^^y'^LT^, such that 

£ D -UL UR + iu^ d). (4) 

Here, we concentrate on possible experimental con¬ 
straints on the diagonal entry Yc = (Y“) 22 . For conve¬ 
nience, we parametrize the deviations from the SM pre¬ 
diction {Cu = Cd = 0 ) in terms of Kq = YqvJ{\f2mq) ^ 
1 , which we assume to be real for simplicityld 

III. THE QCD-YUKAWA pp he PROCESS 

We consider the production of a Higgs boson in asso¬ 
ciation with a charm-quark jet. At the LHC, the main 


^ In the following we assume the top and bottom Yukawa cou¬ 
plings to be constrained close to their SM values after the high- 
luminosity LHC run. 



FIG. 1. Diagrams contributing to pp —>■ he at leading order. 
Black dots correspond to vertices where the Yukawa coupling 
Yc enters, while the crossed vertex corresponds to the SM-like 
top triangle, integrated out. 


partonic process inducing this hnal state is gc —>■ he and 
the corresponding Feynman diagrams are presented in 
Figure The charm Yukawa coupling, depicted as a 
black dot, enters in the first two graphs, that yield a 
contribution to the amplitude of 0{gsYc). The t—channel 
diagram turns out to be largely dominant. The third dia¬ 
gram is formally of higher order in but is enhanced by 
the top-quark Yukawa coupling. Here the crossed vertex 
corresponds to the effective ggh interaction obtained by 
integrating out the top quark. This diagram yields the 
contribution to the amplitude that survives in the limit 
Kc — t 0 (see Table |T|. 

The challenge of the proposed process is to tag the 
charm-quark jet, as in /i —> cc. However, as anticipated, 
it offers some interesting virtues compared to /i —> cc. 
In particular, it allows us to fully reconstruct the Higgs 
boson in a clean decay channel such as > 77 or —>■ 
WW, and it requires only a single charm tag. The main 
drawback is that the process does not vanish in the limit 
Yc —f 0 (contrary to h ^ cc) requiring a good theoretical 
control on the cross section as a function of Yc- While 
a full analysis, including the optimization of the event 
selection, is beyond the scope of this article, here we just 
want to examine the potential of the channel by deriving 
the expected number of signal and background events, 
based on reasonable efficiency assumptions. 

We have calculated the cross section of pp —>■ he at 
leading order in QCD (including the effective ggh as dis¬ 
cussed above) at the LHC with 14TeV center-of-mass 
energy for various values of Key employing MadGraphS 
[lOj . with a tailored model file and CTEQ6L1 parton 
distribution functions. Using uicimz) = 0.63 GeV and 
rrih = 125 GeV, for k^ = I {i.e., the SM) we obtain a cross 
section of cr{pp —>■ he) = 166.1 fb, employing the default 
cuts of ptU) > 20 GeV, g)]) < 5, Ai?(ji, 42 ) > 0.4 for all 
processes considered here. In the following, we focus on 
the —>■ 77 decay channel, with a branching fraction of 
B{h —>■ 77) = 0.0023. This leads to = 2292 events at 
the HL-LHC with 3000 fb“^, taking into account also the 
pp —>■ he process. Assuming a charm-tagging efficiency 
of Cc = 0.4 (see e.g. Ref. M), we finally end up with 
S = CeSo = 917 signal events. The different number of 
events obtained by varying Kc are reported in Table |lj 

The main backgrounds to the process studied here 
are pp —>■ hg, with the gluon mis-identified as a charm 
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2.75 
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3.25 
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3.75 

4 

4.25 
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s 

1097 
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1424 
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TABLE I. Number of Signal events S'(kc) in dependence on 
the charm-quark Yukawa coupling. See text for details. 



Kc 

FIG. 2. The expected p-value for a given value of Kc from 
the process pp —>■ he at the 14 TeV LHC with 3000 fb“^ and a 
conservative assumption for the theoretical uncertainty. See 
text for details. 


quark, as well as pp —>■ hb, with the bottom quark be¬ 
ing mis-tagged. In the first case, we treat separately 
the case pp —>■ hcc, where only one charm-quark jet is 
reconstructed and the case where the gluon produces a 
light quark jet. The backgrounds feature (7{pp —>■ hg) = 
12.25 pb, cr{pp — 7 > hb) = 203 fb, as well as (T{pp —>■ hcc) = 
55 fb. We employ a conservative assumptions for the jet 
reconstruction efficiency of 1 — e™*®® = 95%, as well as 
5 —)■ c and 6 —> c mis-tag rates of eg->.c = 1% and 
eb->.c = 30%. With these figures we obtain B = 1705 
background events at 3000 fb“^, leading to N{kc = 1) = 
S{kc = 1) + B = 2622 total events. We then assume 
a statistical error on the total number of events (VN) 
and a theoretical (relative) error on the signal events of 
20%. The latter is deduced by the recent next-to-leading 
order (NLO) analysis of the Higgs production in associ¬ 
ation with bottom quarks m- Finally, statistical and 
theoretical error are added in quadrature}^ 

In the following, we want to examine the expected 
constraints that can be set on Kc from the process un¬ 
der consideration. To this purpose, we assume the SM 
to be true and calculate how many standard deviations 
AN{kc) away a prediction N{kc) is from N{Kf. = 1 ), 
which is the expected outcome of the experiment. The 
values of Kc that lead to a discrepancy of more than n 
standard deviations are then expected to be excluded at 


® The two dominant backgrounds, pp —^ hb and pp —^ hg, can both 
be directly measured at the LHC with specific tags (inverted b 
vs. c tag for the former and light-quark-jet tag for the latter) - 
this is why we do not assign an additional theory error to them. 


na. We plot the correspondingp-value, p(kc), in Figure}^ 
approximating the Poisson distribution of the number of 
events by a Gaussian. The Ict and 2cr equivalents are 
depicted by the solid and dashed lines, respectively. A 
conservative estimate for the expected l-cr (95% CL) con¬ 
straint on Kc is thus obtained as 

|Kel < 2.5 (3.9), (5) 

which lies in the ballpark of the results quoted in [5], 
where the latter combines ATLAS and CMS to arrive at 
2 X 3000 fb“^ of integrated luminosity. 

On the other hand, an improved prediction of the SM 
cross section a(pp —>■ he), leading to 5th = 10 %, would 
strengthen our expected l-cr (95% CL) limit to 

I«cl < 1-9 (2.6), (6) 

approaching the SM value of Yc- 

We note that optimized cuts can still increase S/B and 
in particular lead to an enhanced sensitivity on Kc- As 
the statistics at 3000 fb“^ is large enough, there are good 
prospects to still improve the bounds. A corresponding 
detailed investigation, including detector simulation, is 
beyond the purpose of this letter and can be performed 
best by the experimental community. 

We further stress that the dominant source of uncer¬ 
tainty, at present, is the theoretical error on a{pp — > he). 
We have indeed checked that the result does not change 
significantly worsening the 5 —>■ c and b ^ c mis-tag rates 
to 5% and 40%, respectively. As far as the reliability (and 
possible reduction) on the theoretical error is concerned, 
a promising possibility would be a dedicated calculation 
of <j{pp —)■ he)/a{pp — )• hb) at NLO (or NNLO), as a 
function of YcjYb, supplemented by measurements of this 
ratio and (T{pp —> hb) with a combination of normal and 
inverted b vs. c tags. 

IV. THE ELECTROWEAK pp hM PROCESS 

As anticipated in the introduction, the production of 
the Higgs boson in association with charm can proceed 
also via electroweak interactions, starting form an initial 
charm-less qq' state {ud —>■ hW^*'> —>■ hes). The case of 
an on-shell W producing a charm jet can be discrimi¬ 
nated from the QCD-Yukawa process by means of ap¬ 
propriate cuts on the jet momentum. Less obvious is the 
discrimination in the case of a virtual W* producing a 
low-momentum c-jet, or even a single charmed hadron. 
In the following we estimate in detail the specific case of 
the single meson production: pp hM., with M being 
a charmed meson or a charmonium state. 

The leading partonic amplitude within the SM is 
shown in Fig. Following Refs. mill], we parameterize 
the quark currents appearing in the initial and final state 
with arbitrary vector and axial couplings: 

Jq,^j = Q\ 9 v,^J 1 ^ + 9 A,ij 7 ^ 75 ) 9 ^ ■ 


(7) 
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FIG. 3. Diagram contributing to pp —>■ hM at leading order, 
where V = W^,Z. 


The matrix element of the current that generates the me¬ 
son in the final state assumes one of the following struc¬ 
tures, depending on the spin of Ai: 


{M{p,e)\J:;\0) 


1 

2 


9 v fv p'^ 

9v fv mp £'' 


M = V 
M = V, 


( 8 ) 


where Jm is the meson’s decay constant, and pM encodes 
the dependence on the coupling to the relevant gauge 
boson {g-p = gA,ij, gv = gv,ij for a {q^g'') meson) With 
this notation, the SM expression for the partonic cross 
section for the case of a pseudoscalar meson reads 


a{qq' -A hV)sM{q^) 


b7(S'Kv‘^{q^ — myY ’ 


(9) 


where V = , Z, and we have suppressed the indices of 

9 a,V for simplicity. The vector case has the same func¬ 
tional form with 7^ —>■ V, up to tiny O (my/wy) cor¬ 
rections. In the above expression, q"^ denotes the total 
momentum of the initial state in the partonic process 
and 


A(a = Wl -2 


mt 


m 




Convoluting the cross sections with the appropriate PDF 
in the region 130 < \/^ < 1 TeV, and assuming an in- 
tegrated luminosity of 3000 fb“^ we obtain the expected 
number of events for each channel at the HL-LHC. The 
results, summarized in Table |Tlj show that these pro¬ 
cesses will not be observable at the SM level, and cer¬ 
tainly do not represent a dangerous background for the 
QCD-Yukawa process discussed in Sect. |III| 

Given the smallness of the SM signal, it is worth to 
investigate if these cross sections can be significantly 
altered beyond the SM. This can be done generalizing 
the approach of Refs. m m- The leading (helicity- 
conserving) transition amplitude can be decomposed in 
full generality as 


Aiq^f ^ hM) = -J^T^,{M\J^\0) . (11) 


gv,uu — (4 3 ^w)’ 9 v,dd — ( 4 3 *w) > 9A,uu — 

9A,dd 4c^ ’ 3v,u^dd 9A,u^d^ 2\/2 


Ghannel 

"m-M /a4 

(MeV) 

Events @ HL-LHG 

Method (a) Method (b) 

9c 

2984 200 

0.10 

0.08 

JM 

3100 410 

0.08 

0.07 

Dt 

1968 250 

0.48 

0.40 

d:± 

2112 325 

0.84 

0.69 


TABLE II. Expected number of hM associated production 
events at HL-LHC (14 TeV and 3000 fb~^) in the energy 
region 130 < ^/q^ < 1 TeV for representative charmed- 
meson final states. The results reported under Method (a) 
are obtained rescaling bin-by-bin the cross section distribu¬ 
tion of Drell-Yan processes provided by MadGraph 5 m- 
The computation of Method (b) is performed via numerical 
convolution of the analytic cross section with the PDF of the 
MSTW 2008 libraries [13]. Both account only for SM contri¬ 
butions. 


To a good accuracy the quark current is conserved 
[q^iJq = 0), and the tensor can be decomposed in 
terms of only four Lorentz structures. Using the same 
notation as in Ref. [5]: 

= fi{q^)9t^v + f2{q^)Pti.Pv + /3(g^)(p • 7S/xi' - PMu) 
+fi{q^)£tiupaP^q'" , ( 12 ) 

where q^ is the total momentum of the quark pair in the 
initial state, and p^ is the meson momentum (p^ = to^). 
With these notations the partonic cross section reads 

a{qq' hr){q^) = {py + 9 a) 

X |/i('7^) + w^/ 2 (g^)|^ q^>?{q^) , (13) 

where, similarly to the SM case, (j{qq' -A hV) has the 
same functional form up to tiny O {rriy/rn^) corrections. 
Neglecting the latter terms, we obtain 

a{qq' -A hM)BSM , 2 \ 
a{qq' -A hM)sM ^ 

where ff^{q^) oc l/(v{q^ — my)) and we disregard poten¬ 
tial changes to the fermionic currents. Deviation from the 
SM are thus induced by possible non-pole-terms (i.e. con¬ 
tact terms) in the form factor fi{q^). Within a generic 
effective-field theory (EFT) approach to Higgs physics 
(both linear and non-linear EFT), contact terms in fi{q^) 
are generated by dimension-six operators. However, their 
effect would show-up exactly in the same functional form 
either in the on-shell associated production {pp -A Vh) 
or in ft, —>■ VM decays, that share the same current struc¬ 
ture [Sill]. Since the latter processes can be measured 
(or at least bounded) to a better accuracy, we conclude 
that <j{pp -A hM) is not a very sensitive probe of generic 
extensions of the SM. 


h{9^) 


hm 


SM 


(14) 
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V. CONCLUSIONS 


In this letter, we proposed a new strategy for the mea¬ 
surement of the Yukawa coupling of the charm quark: the 
measurement of the production cross section of the Higgs 
boson in association with a charm jet. A first estimate 
showed that could be determined at a level approach¬ 
ing the SM value in this channel, which offers virtues 
and drawbacks quite different with respect to the > cc 
search. A fully realistic analysis was beyond the scope of 
the present paper. A more realistic evaluation of the effi¬ 
ciencies is likely to decrease the number of signal events 
S compared to our naive estimate; however, as we have 
discussed, the sensitivity on could even increase with 
properly designed b and c tag strategies aimed to measure 


the background from data and to reduce the theoretical 
error on the normalization of the cross-section. This first 
analysis therefore calls for more detailed studies both on 
the theory and on the experimental side. 
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